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In the presence of a linear potential with an arbitrary time-dependence, Hirota method is de- 
veloped carefully for applying into the effective mean-field model of quasi-one-dimensional Bose- 
Einstein condensation with repulsive interaction. We obtain the exact nonautonomous soliton so- 
lution (NSS) analytically. These solutions show that the time-dependent potential can affect the 
velocity of NSS. In some special cases the velocity has the character of both increase and oscillation 
with time. A detail analysis for the asymptotic behaviour of solutions shows that the collision of 
two NSSs is elastic. 
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I. INTRODUCTION 



The concept of soliton was introduced firstly by 
Zabusky and Kruskal j| to characterize nonlinear soli- 
tary waves that do not disperse and completely preserve 
their localized form and speeds during propagation and 
after a collision. This intrinsic favorable property of soli- 
ton has motived a great of attention on the nonlinear 
systems in many fields of physics, especially in high-rate 
telecommunications with optical fibers and condensate 
physics. Hasegawa and Tappert 2] derived the nonlin- 
ear Schrodinger (NLS) equation model in fiber and firstly 
predicted optical soliton, and then experimental verifica- 
tion has successfully been carried out by Mollenauer et 
al. Q. Since then, optical solitons have been the objects 
of extensive theoretical and experimental studies in the 
past three decades for their potential applications in long 
distance communication and all-optical ultrafast switch- 
ingdevices. Recently, the controllable soliton solutions 
are of interest in the field of nonlinear optics and 
condensate physics, and then the term of nonautonomous 
solitons [6] introduced firstly. In fact, different aspects of 
dynamics in nonautonomous models J8i,i9|| in linear poten- 
tials have been investigated theoretically. Strictly speak- 
ing, this nonautonomous solutions obtained could not be 
considered as canonical solitons. Fortunately, the real- 
ization of the Bose-Einstein condensation (BEG) [13, Ell 
offered a good examples of the nonautonomous systems 
in condensate physics. 

With the realization of BEGs the exploration of the 
nonlinear properties of matter waves has been paid more 
particular interest. One of them is macroscopic ally ex- 
cited BEGs, such as vortices 12| and solitons jl3l - [l9| . At 
zero temperature the dynamics of BEG is well described 
by the time-dependent Gross-Pitaevskii (G-P) equation, 
and the nonlinearity results from the interatomic interac- 
tions. Depending on the attractive or repulsive nature of 
the interatomic interactions, G-P equation has of cither 
bright or dark soliton solutions, respectively. A bright 
soliton [20l - [2^ in BEG is expected for the balance be- 



tween the dispersion and the attractive mean-field en- 
ergy. However, large condensates are necessarily associ- 
ated with repulsive interaction, for which bright soliton 
might seem impossible because the nonlinearity cannot 
compensate for the kinetic energy part in the atomic dy- 
namics. So it is interesting to explore the property for 
dark soliton of BEG. A dark soliton [H [l!] in BEG is 
a macroscopic excitation characterized by a local density 
minimum and a phase gradient of the wave function at 
the position of the minimum. Under the different con- 
ditions many soliton solutions [15nl9l HE HI] have been 
obtained, as well as the dynamics of the excitation of the 
condensate was discussed. When the longitudinal dimen- 
sion of BEG is much longer than its transverse dimension 
which is the order of its healing length, the G-P equation 
can be reduced to the quasi-one-dimensiona l (q uasi-lD) 
regime. This trapped quasi-one-dimensional [27[ conden- 
sate has offered an useful tool to investigate the nonlin- 
ear excitations such as solitons and vortices, which are 
more stable than in 3D, where the solitons suffer from the 
transverse instability and the vortices can bend. So the 
studies of both theory and experiment are very impor- 
tant for the soliton excitation in quasi-one-dimensional 
BEG. 

The effective mean-field model of a quasi- ID BEG in 
a linear potential with an arbitrary time-dependence is 
given by 



in 'I' = 

dT 2m dX^ 



■^ + Xf{T)-^ + g\^\'-^, (1) 



where 



I^E'I dX — N \s the number of atoms in the 



condensate. The interacting constant of two-atom '\s g = 
2h?a/mlf^ [15, where m is the mass of the atom, a is the 
s-wave scattering length (a < for attractive interaction; 
while a > for repulsive interaction), and Zo = \/fi/mu]Q 
is the characteristic extension length of the ground state 
wave function of harmonic oscillator. For pithiness, we 
introduce x = X/Iq, t = T/ml'^/h, and ijj — ^'/\/iV7o, 



2 



and then Eq. ([T]) reduces to the dimensionless form 



2d^ 



7/- + a;/(t)V^-Ml^N = 0, (2) 



2Nloa, and f {t) 



ml} 



/( 



T ■ 



As in 



where ^ - ^i»iuu,, miu j v--; - — w 
Ref. [6], Eq. ([2]) was called nonautonomous NLS model 
in linear potential. When / (t) = constant and a < 0, 
the Lax pair and NSSs in inhomogeneous plasma has 
been constructed by the inverse scattering method 
F-expansion method [2^ and Hirota method [3(i] were 
also developed to construct the bright NSSs in quasi-lD 
BECs. However, the dynamic property of NSS of Eq. 
hasn't well explored, and Hirota method developed 
for Eq. ([2]) with the repulsive interaction is also very 
interesting. 

In the present paper we consider mainly the dynamics 
of dark nonautonomous soliton which can be affected by 
adjusting the external linear time-dependent potential. 
In the following section we demonstrate how to construct 
the exact dark NSSs of Eq. and the corresponding 
properties of such solutions are studied in detail. 



II. DEVELOPED HIROTA METHOD AND ONE 
NONAUTONOMOUS SOLITON SOLUTION 

Hirota method [3l| is an effective straightforward tech- 
nique to solve the nonlinear equations. In order to clear 
the derivation of solution we introduce the main idea of 
Hirota method briefly. Firstly, it apply a direct trans- 
formation to the nonlinear equation. Then, by means 
of some skillful bilinear operators the nonlinear equation 
can be decoupled into a series of equations. With some 
reasonable assumptions the exact solutions can be con- 
structed effectively. However, in the presence of the time- 
dependent potential the application of Hirota method 
should be more careful to get NSSs of Eq. (0) in the 
case of the repulsive interaction. 

Performing the normal procedure, we consider the 
complex function G (x, t) and the real function F {x, t) 
forming the transformation 



Gix,t) 



F{x,t)- 

Substituting Eq. ^ into Eq. Q we have 



(3) 



FiiDt + ^)G-F + GF^xf{t)-G{^F-F + ^j,GG) = 0, 

(4) 

where the overbar denotes the complex conjugate, Dt and 
are called Hirota bilinear operators defined by 



D^^D^G{x,t)-F{x',t') 



(5) 



In the absence of the time-dependent potential, i.e., 
^f{t) = 0, Eq. dl]) can be decoupled easily into two 
equations. In order to get exact dark NSSs of Eq. ([2]), 
a real parameter A to be determined should be added 
to Hirota bilinear operators in the presence of the time- 
dependent potential, and then many attempts show that 
Eq. dl]) can be decoupled into 



AiG-F = 0, A2F ■ F = -nGG, 



(6) 



where the overbar denotes the complex conjugate, and 
Hirota bilinear operators Ai and A2 are given by 

Ai=iDt + ^Dl + xf{t)- X, 

M = \dI - A. 

The derivation of Eq. ^ has made Eq. ^ into the 
normal procedure of Hirota method. The spatial and 
time dependence term xf{t) will play an important role 
for getting the exact NSSs as shown later. 

If the expressions of G {x, t) and F {x, t) are obtained 
from Eq. the exact dark NSSs can be expressed 

analytically. For this purpose we assume that 



G = Go(l + xGi),F=(l + xFi) 



(7) 



where x is an arbitrary auxiliary parameter which will 
be absorbed in the expression of NSSs. Substituting Eq. 
([7]) into Eq. ^ , and collecting the coefficients with same 
power of X, we have 

(1) for the coefficient of x° 



0, 



(8) 
(9) 

(10) 
(11) 



Ai (Go • 1 
i2(l-l) = -/iGoGo, 

(2) for the coefficient of 

Ai (GoGi-l-f Go-Fi) = 0, 

A2{Fi-l + l- Fi) = -fiGoGo (Gi + Gi 

(3) for the coefficient of 

Ai (GoGi • Fi) =. 0, (12) 
A2{Fi-Fi)^-fiGoGoGiGi. (13) 

Using the definition of Hirota bilinear operator ([S]) the 
above equations can be expressed in detail. Considering 
the presence of the term xf{t) in Eq. ^ we assume Gg 
has the form 



where 



Go = 70 e*' 



Tjo = Po{t)x + no{t), 



(14) 



(15) 



with Po(t) and ^o(t) is to be determined, respectively. 
Substituting Go into Eq. ([5]) we have 



x—x'^t—t' 



= [-PoAt) + fit)] ^ - ^oAt) ~ 0^0 (^) - 
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which imphes the solution 



Poit) 



f{T)dT + ^0, 



1 /■* 



(16) 



where and (^o is an arbitrary real constant, respectively. 
Prom the restriction of Eq. ([9]) we get |7o|^ = A//i which 
shows that the exist of dark NSS demand the parameter 
A > 0. For convenience 70 can be chosen as 70 = \f\f\i. 

Expanding Eqs. (|TU| and (TTlT) with the definition of 
Eq. ([5]) one can see that G\ and F\ admit the expression 



Gi = Zx exp 771, Fi = exp 771, 



(17) 



where the parameter Z\ to be determined is complex, 
and the real parameter r]\ is given by 

771 = Pi (t) X + (t) . (18) 

Substituting Eqs. (O and ((l2l) into Eq. we have 

= z (Zi - 1) (Pi,tx + rii^t + PoA) + ^ (^1 + 1)-Pi'. 

(19) 

In the case of Z\ ^ ±l(the case oi Z\ = ±1 will be 
discussed in below), the solution of Eq. ([T9|) is P\^t = 0, 
i.e.. Pi should be independent on and Zi is given by 



Zi 



iPaPi - ip2 + ifli. 



iP^Pl + ip2 + lili^t ' 



(20) 



which shows that \Zi\ 
get the restriction 



1. From Eqs. ((Ill) to (US]) we 

(21) 



Pf = 2A - A (Zi + Zi) . 
Then from Eq. and Eq. dUl) we obtain 



^1 



V4A - P^ + iPi 
V4A - P2 - iPi ^ 



r!i = -Pi ^ Po (r) dT + ^^iX-P^Pit + Ci, 



(22) 



where Ci is a real constant. Pi is an arbitrary real pa- 
rameter, and A > Pi'^/4. 

With the help of Eqs. ^ and (d?]), after ab- 

sorbing X, the exact dark NSS of Eq. ([2]) can be derived 
as 




(1 + Zi)- (l-Zi)tanh 



(23) 



where the parameters 770, 771, and Zi have been given in 
Eqs. (dil), (HH), (UHl) and 

If one chose A = Ar/^ and Pi = 4a 17;^, the solution in 
Eq. ((23)) has the same form as the general solution ob- 
tained in the framework '9]. In the case of f{t) = 0, the 



solution in Eq. (|23p reduces to the one soliton solution 
of the normal NLS equation. When f{t) = constant, the 
solution ((23]) represents the nonlinear dark wave propaga- 
tion in linearly inhomogeneous plasma or optical fibre 
with the abnormal dispersion. As f{t) — bi + I cos (utt), 
the solution ((25)) denotes NSS in BEG with considering 
the coupling of the external field and the effect of gravity. 

From the NSS in Eq. we clear two special case 

mentioned before, i.e., Zi = ±1. When Zi = 1, the s olu- 
tion ((25)) reduces to plane- wave solution -ipi ~ y^X/fie'^° , 
which corresponds to the uniform distribution density of 
bosons. On the other hand, when Zi = — 1 the solu- 
tion (125)) becomes ipi = —y/X/fj,e^^° tanh (7/1/2), where 

771 — 2y/Xx — 2\/X J* Po (r) dr + (i. This solution repre- 
sents the black NSS of BEC in a linear potential with an 
arbitrary time-dependence which is caused by the cou- 
pling of external field and the effect of gravity. 

We also find the effect of term xf{t) from the NSS in 
Eq. (|23p . As shown in the expression of 770, the term 
xf{t) can only contribute a phase to the background. 
The width of nonautonomous soliton, defined by 1/Pi, 
is not affected by the time-dependent external potential. 
From Eqs. ((T8)) and ((22)) we get the velocity 



which shows that the time-dependent potential play 
an important role for the velocity of dark NSS. When 
f{t) = bi + lcos{ujt), the velocity becomes Vi = bit + 



I /uj sin (ujt) + Co ^ Y A — (Pi/2) , which increases and os- 
cillates with time. These results shows that under the 
effect of gravity, the BEC slides down where the dark 
nonautonomous soliton slides down and oscillates in time. 
It can be realized by controlling the external field. 



III. COLLISION OF TWO DARK 
NONAUTONOMOUS SOLITONS 

In this section we will give the analytical expression 
of two dark NSSs of Eq. ((2]). In this case G{x,t) and 
F {x, t) of Eq. ((3)) are assumed as 

G = Go (1 + xGi + x'G2) , P - 1 + x^^i + (24) 

where Go has been obtained in Eq. ()14p . Employing 
the similar procedure of the above section we obtain the 
following set of equations from Eq. ((Bj) 

(1) for the coefficient of x 

£i(Gi-l + l-Pi) = 0, (25) 
i2(l-Pi+Pi-l)--A(Gi + Gi), (26) 

(2) for the coefficient of 

£i (1 • P2 + Gi • Pi + G2 • 1) = 0, (27) 

A2 (1 • P2 + Pi • Pi + P2 • 1) = -A (G2 + GiGi + G2) , 

(28) 
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(3) for the coefficient of 

£i{Gi-F2+G2-Fi) = 0, (29) 

A2 (Fi ■F2+F2- Fi) = -A (G2G1 + G1G2) , (30) 

(4) for the coefficient of 

£iG2 • F2 - 0, (31) 

i2^^2 .F2 = -AG2G2, (32) 

where £1 — iDt + ^D-^+Po{t)Dx, and A2 is given before. 



where the real parameter A12 is given by 



It is obvious that one can solve the equations ((25|) to 
(15^ in turn with the reasonable expressions of Gi and 
Fi. A detail analysis shows that Gi and Fi admit the 
forms 

Gi = Zi exp 771 + Z2 exp 772 , -Fi = exp rji + exp 772 , (33) 

where Zj is complex and rij has the form 

rjj^P^(t)x + n,{t) J ^1,2, (34) 

with the parameter Pj(t) and f2j(t) is to be determined, 
respectively. Substituting Eq. into Eq. (^51) we have 

= e^Ml^i - mPi.tx + ^PoPi + i^i,t) + ^^Y^Pi] 

+ €"^[(^2 - l){lP2,tX + iPoF2 + ^^^2,^) + 

In the case oi Zi ^ ±.1 and Z2 7^ ±1, the above equation 
implies that Pj t {t) = 0, j = 1,2, i.e., Pj is independent 
on t, and Zj is given by 



(35) 



which shows that \Zj\ = 1, j — 1,2. From Eq. (^5)) we 
have the restriction 

F/ = 2A-A(Z,+Z,),j-l,2. (36) 

From Eqs. ((35)) and (p6| we have 



4A - P/ + ,p^. 
4A - P| - iPj 



-P, ^* Po (r) dr + i y^4A - PfP.t + Q , (37) 



where Q is a real constant, P, is an arbitrary real pa- 
rameter, and A > F//4, j — 1, 2. 

Substituting Eq. §3^ into Eqs. ([221) and after 
a tedious and expatiatory calculation we obtain the ex- 
pressions of G2 and P2 as 

G2 = A12Z1Z2 exp (771 + 772) , P2 = A12 exp (771 -I- 772) , 

(38) 



4A - P1P2 - v/4A - PiV4A - P| 

^12 — — 

4A + P1P2 - V4A-PiV4A-p2 



Now we have obtained the expression of Go,Gi,G2,Pi 
and P2 in Eq. With the help of Eqs. (155)) to 

one can find the Eqs. (|25)) to (|5^ are satisfied to the 
moment. 

With Eqs. Q, (O, ([33]), and dSl]), while absorbing 
the parameter x, we obtain the dark two NSS of Eq. ^ 
as 



^2 



/ A^,,^^ 1 + Zie^i + Z2e''^ + ^12^1^26"^+"^ 



(39) 

When f{t) ~ 0, the solution (15^ denotes dark two soli- 
tons interaction of the normal NLS equation. When 
f{t) =constant, the solution (I39p represents the dynam- 
ics of two nonautonomous nonlinear waves in linearly in- 
homogeneous plasma or optical fibre with the abnormal 
dispersion. As shown before, the expressions (1231) and 
([M| imply that Hirota method has more advantage for 
getting such solutions as well. 

The solution in Eq. (|39)) describes a general scattering 
process of two dark NSSs with different center velocity 
Vi and V2, respectively. From Eqs. (|M)) and (|57)) we get 
each velocity as 



/' 

Jo 



V, = / /(r)dr-feo 



4A-P2,j = 1,2. 



Under the proper parameters two NSSs can move toward 
each other, one with the velocity Vi, while the other with 
V2. In order to understand the nature of two NSSs inter- 
action, we analyze the asymptotic behave of the solution 
in Eq. (|39)) . Asymptotically, the solution in Eq. (|39)) can 
be written as a combination of two NSSs in Eq. (|23l) . 
The asymptotic form of two NSSs in limits t —00 and 
i — )■ 00 is similar to that of the one NSS in Eq. (05)) . 

(i) Before collision (limit t ^ -co) 

(a) Nonautonomous soliton 1 (771 w 0, 772 — > —00) 



;H-Zi)-(1-Zi)tanh|- , (40) 




(b) Nonautonomous soliton 2(7;2~0, 7/1— >oo) 



(1 + Z2) - (1 - Z2) tanh - (772 + <5o) 




(ii) After collision (limit t 00) 

(a) Nonautonomous soliton 1 (771 « 0, 772 — 00) 



1 



(41) 




(1 + Zi) - (1-Zi) tanh- (7/1 +(5o) 

(42) 
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(b) Nonautonomous soliton 2 ( 772 ~ 0, 771 — ^ —00) 




(1 + Z2)-(1-Z2)tanh^ , (43) 



where the center shift is given by Sq — Inv4i2. By analyz- 
ing the asymptotic behave in detail we know that there 
is no change of the amplitude for each NSS during col- 
lision, while one should notice that the factor \Zj\ — 1, 
j = 1,2, again. However, from Eq. (gOl) to Eq. 
we find a phase exchange 60 for soliton 1 and soliton 2 
during collision. These results show that the collision of 
two NSSs is elastic. 



trary time-dependence, and Hirota method is also devel- 
oped. With the skillful assumption the exact dark NSSs 
are constructed effectively. From these results we find 
the time-dependent potential can affect the velocity of 
NSS. In some special cases the velocity of NSS in quasi- 
one-dimensional BEC increases and oscillations in time, 
the BEC slides down under the effect of gravity. We also 
investigate the asymptotic behave of two NSSs which de- 
notes the elastic collision. 



ACKNOWLEDGEMENT 



IV. CONCLUSION 

In this paper, we report the exact dark NSSs of quasi- 
one-dimensional BEC in a linear potential with an arbi- 



This work is supported by NSF of China under Grant 
No. 10874038, the Natural Science Foundation of Hebei 
Province of China under Grant No. A2008000006, and 
the key subject construction project of Hebei Provincial 
University of China. 



[1] N. J. Zabusky, and M. D. Kruskal, Phys. Rev. Lett. 15 
(1965) 240. 

[2] A. Hasegawa, F. Tappert, Appl. Phys. Lett. 23 (1973) 
142. 

[3] L.F. Mollenauer, R.H. Stolen, J.P. Gordon, Phys. Rev. 

Lett. 45 (1980) 1095. 
[4] V. N. Serkin and A. Hasegawa, Phys. Rev. Lett. 85 (2000) 

4502. 

[5] V. N. Serkin, et al., JETP Lett. 74 (2000) 573; 
V. N. Serkin, JETP Lett. 72 (2000) 89; 
V. N. Serkin, et al, IEEE J. Quantum Electron. 8 (2002) 
418. 

[6] V. N. Serkin, A. Hasegawa, and T. L. Belyaeva, Phys. 

Rev. Lett. 98 (2007) 074102. 
[7] Lei Wu, Jie-Fang Zhang, Lu Li, C. Finot, and K. 

Porsezian, Phys. Rev. A 78 (2008) 053807 
[8] H. H. Chen and C. S. Liu, Phys. Rev. Lett. 37 (1976) 

693; Phys. Fluids 21 (1978) 377 
[9] V. V. Konotop, Phys. Rev. E 47 (1993) 1423; 

V. V. Konotop, O. A. Chubykalo, and L. Vazquez, Phys. 

Rev. E 48 (1993) 563. 

V. V. Konotop, Theoretical and Mathematical Physics, 
99 (1994) 687. 

[10] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 
Wieman, and E. A. Cornell, Science 269 (1995) 198; 
K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van 
Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, 
Phys. Rev. Lett. 75 (1995) 3969. 

[11] C. C. Bradley, C. A. Sackett, J. J. ToUett, and R. G. 
Hulet, Phys. Rev. Lett. 75 (1995) 1687; 78 (1997) 985; 
D. G. Fried, T. C. KiUian, L. Willmann, D. Landhuis, 
S. C. Moss, D. Kleppner, and T. J. Greytak, Phys. Rev. 
Lett. 81 (1998) 3811. 

[12] R. J. Donnely, Quantized Vortices in Helium II (Cam- 
bridge University Press, Cambridge, 1991). 
M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 
Hall, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 
83 (1999) 2498. 



A. L. Fetter and A. A. Svidzinsky, J. Phys.: Condens. 
Matter 13 (2001) R135. 
[13] Y. S. Kivshar and B. Luther-Davies, Phys. Rep. 298 
(1998) 81. 

[14] P. A. Ruprecht, M. J. Holland, K. Burnett, and M. Ed- 
wards, Phys. Rev. A 51 (1995) 4704. 

[15] W. Zhang, D. F. Walls, and B. C. Sanders, Phys. Rev. 
Lett. 72 (1994) 60; 

W. P. Reinhardt and C. W. Clark, J. Phys. B 30 (1997) 
L785; 

A. D. Jackson, G. M. Kavoulakis, and C. J. Pethick, 

Phys. Rev. A 58 (1998) 2417. 
[16] A. E. Muryshev, H. B. van Linden van den Heuvell, and 

GV. Shlyapnikov, Phys. Rev. A 60 (1999) R2665; 

Lincoln D Carr, Mary Ann Leung, and William P Rein- 
hardt, J. Phys. B 33 (2000) 3983. 
[17] J. Denschlag, et al.. Science 287 (2000) 97. 
[18] J. C. Bronski, L. D. Carr, B. Deconinck, and J. N. Kutz, 

Phys. Rev. Lett. 86 (2001) 1402. 
[19] L. D. Carr, C. W. Clark, and W. P. Reinhardt, Phys. 

Rev. A 62 (2000) 063610; 62 (2000) 063611; 

D. L. Feder et al., Phys. Rev. A 62 (2000) 053606. 
[20] Kevin E. Strecker, Guthrie B. Partridge, Andrew G. 

Truscott, and Randall G. Hulet, Nature (London) 417 

(2002) 150. 

[21] L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. Cu- 

bizolles, L. D. Carr, Y. Castin, and C. Salomon, Science 

296 (2002) 1290. 
[22] U. Al Khawaja, H. T. C. Stoof, R. G. Hulet, K. E. 

Strecker, and G. B. Partridge, Phys. Rev. Lett. 89 (2002) 

200404. 

[23] R. Dum, J. I. Cirac, M. Lewenstein, and P. ZoUer, Phys. 

Rev. Lett. 80 (1998) 2972. 
[24] S. Burger, K. Bongs, S. Dettmer, W. Ertmer, and K. 

Sengstock, Phys. Rev. Lett. 83 (1999) 5198. 
[25] Shun- Jin Wang, Cheng-Long Jia, Dun Zhao, Hong-Gang 

Luo, and Jun-Hong An, Phys. Rev. A 68 (2003) 015601. 
[26] Zai-Dong Li, P. B. He, Lu Li, J.-Q. Liang, and W. M. 



6 



Liu, Phys. Rev. A 71 (2005) 053611; 
Lu Li, Zaidong Li, Boris A. Malomcd, Dumitru Miha- 
lache, and W. M. Liu, Phys. Rev. A 72 (2005) 033611; 
Q. Y. Li, Z. W. Xie, L. Li, Z. D. Li, and J. Q. Liang, 
Annals of Physics 312 (2004) 128. 
[27] A. Gorlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, T. 
L. Gustavson, and J. R. Abo-Shaeer, Phys. Rev. Lett. 87 
(2001) 130402; 

S. Dettmer et al., ibid. 87 (2001) 160406. 



[28] D. S. Petrov, G.V. Shlyapnikov, and J. T. M. Walraven, 

Phys. Rev. Lett. 85 (2000) 3745. 
[29] Qin Yang, and Jie-fang Zhang, Optics Communications 

258 (2006) 35. 

[30] Zai-Dong Li, Qiu-Yan Li, Xing-Hua Hu, Zhong-Xi Zheng, 

Yubao Sun, Ann. Phys. (N.Y.) 322 (2007) 2545. 
[31] R. Hirota, J. Phys. Soc. Jpn. 51 (1982) 323. 



